arXiv: 1507.04827v3 [cond-mat.str-el] 19 Aug 2015 
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The single-particle spectral properties of the two-dimensional t-J model in the parameter regime 
relevant to cuprate high-temperature superconductors are investigated using cluster perturbation 
theory. Various anomalous features observed in cuprate high-temperature superconductors are 
collectively explained in terms of the dominant modes near the Mott transition in this model. 
Although the behavior of the dominant modes in the low-energy regime is similar to that in the 
two-dimensional Hubbard model, significant differences appear near the Mott transition for the 
high-energy electron removal excitations which can be considered to primarily originate from holon 
modes in one dimension. The overall spectral features are confirmed to remain almost unchanged as 
the cluster size is increased from 4 x 4 to 6 x 6 sites by using a combined method of the non-Abelian 
dynamical density-matrix renormalization group method and cluster perturbation theory. 

PACS numbers: 71.30.-|-h, 71.10.Fd, 74.72.Kf, 79.60.-i 


I. INTRODUCTION 


Cuprate high-temperature (high-Tc) superconductors, 
which are obtained by doping Mott insulators containing 
Cu 02 planes^ exhibit various features that appear 
anomalous from conventional viewpoints<^4i^ Because 
the anomalous features are considered to be related 
to high-Tc superconductivity, the effects of electronic 
correlations near the Mott transition in two-dimensional 
(2D) systems have attracted considerable attention. 

In particular, through an analysis of electronic 
correlations among relevant Cu and O orbitals, the 2D t- 
J model has been derived as a minimal model of high-Tc 
cupratesii^ However, its properties and their relationship 
to the anomalous features are not well understood, 
primarily because of difficulties in dealing with the no- 
double-occupancy constraint. 

Although this model can also be derived effectively 
from the 2D Hubbard model in the large-repulsion 
limit^ it is not clear how similar the two models’ 
properties are in the parameter regime relevant to high- 
Tc cuprates. In fact, the 2D t-J and Hubbard models 
have frequently been studied from different viewpoints: 
the former has been considered a doped Mott insulator 
accessible from slave-particle theorieS j^^d^ while the 
latter has been considered a strongly interacting electron 
system accessible from Fermi liquid theory.^ In some 
studies, double occupancy, which is excluded in the t-J 
model, has been considered important to the anomalous 

features 1^421 

In this paper, by using cluster perturbation theory 
(CPT)^i2£ similarities and dissimilarities in the spectral 
features of these models are clarified. In addition, various 
anomalous features observed in high-Tc cuprates are 
collectively explained in the 2D t-J model, which is an 
effective model of the Cu 02 plane and has no double 
occupancy. A method to reduce cluster-size effects is 
also introduced. 


II. MODELS AND METHODS 
A. Models 

We consider the 2D t-J model defined by the following 
Hamiltonian for t > 0 and J > 0: 

n = -tY^ (cJ,,Cj-„-FH.c.)-FJ^(S'rS'j-in^nj)-Ai^n„ 

{i,j) i 

where Ci^a = (1 - rii-a)ci^a and for the 

annihilation operator Ci^a and number operator of 
an electron with spin cr at site i. Here, Si denotes the 
spin operator at site i, and (i, j) means that sites i and j 
are nearest neighbors on a square lattice. At half-filling 
(doping concentration 6 = 0), the model reduces to the 
Heisenberg model. The t-J model can also be effectively 
obtained for J = AP/U by neglecting the three-site 
ter m — in the laige-U/t limit of the Hubbard model 
defined by the following Hamiltonian: 

"Hnub =-t X 

t i 

In the ground state near the Mott transition, 
ferromagnetic fluctuations might be dominant in the very 
small-J/t regim 022 r.S£ and phase separation occurs in 
the large-J/t regime : Here, focusing attention 
on the parameter regime relevant to high-Tc cuprates 
{J/t K, 0 5) ^13,16 |;jjg CPT results exhibit no 

indication of phase separation [Fig. [Il[m)], we study the 
spectral function defined as A{k,uj) = —ImG(fc,a;)/7r. 
Here, G{k, oj) denotes the retarded single-particle Green 
function for momentum k and frequency co at zero 

temperature— iiid 

B. Methods 

In this paper, CPT is employed: G{k,uj) is 
calculated by connecting cluster Green functions 
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FIG. 1: Results for the 2D t-J model obtained using (4 x 4)-site CPT [(a)-(d) and (f)-(p)] and (6 x 6)-site CPT [(e); yellow 
symbols in (i), (m), and (n) for J/t = 0.5 at 5 = 1/18). (a)-(d) A{k,uj)t for J/t = 0.5 at <5 = 0 {n/t = 4), 0.03, 1/18(« 0.06), 
and 0.2 (from above). The rightmost panels show p{uj)t[= J dkA(k,uj)t/(2Tv)'^] in d(= 2) dimensions. The dotted green line 
in (a) indicates ui at the top of the band. The solid green lines in (b)-(d) indicate u = 0. (e) The same as the left panels of 
(c) but using (6 x 6)-site CPT. The blue numbers indicate dominant modes, (f)-(h) n{k)[= ^du)A[k,uj)] for (b)-(d) for 

a) < 0 [n_(fc)] (solid blue curves with hatches), oi > 0 (dotted red curves), and all uj [nt(fc)] (solid brown lines), (i) Spectral 
weight W[= ^dujp(uj)] for oi < 0 [W-] (red circles), cj > 0 [W+] (blue squares), and all u (green crosses) at J/t — 0.5. 

(j)-(l) The same as in (b)-(d) but for w ~ 0. (m) Chemical potential p, for J/t = 0.5, 0.4, and 0.3 (from above), (n) Energy 
of mode i' at (7r,7r) [e(7r,7r)] (blue squares) for J/t — 0.5. (o) £(o,o) (purple squares) and £(.,r, 7 r) (yellow triangles). The solid 
green line indicates J2vs/t {vs ~ l.lSv^J [H). (p) £p (purple squares) and £(.,r,o) (yellow triangles). The solid green line is 
proportional to J. The red circles in (n) and blue diamonds and red circles in (o) and (p) indicate corresponding energies in 
the 2D Hubbard model {U/t = ‘it/J), taken from Ref. [^ . Gaussian broadening with a standard deviation of O.lt is used. 


through the first-order hopping process i^^'^^i'^°~ — By 
considering superclusters, G{k,uj) at arbitrary <5 can 
be obtained.— In the large-cluster limit, CPT becomes 
exact. Here, to calculate the (4 x 4)-site and (6 x 
6)-site cluster Green functions, exact diagonalization 
and the dynamical density-matrix renormalization group 
(DDMRG) method^ are used, respectively. In the 
DDMRG calculation, the 11(1)0811(2) basis^i^ is 
employed^ and 1000 density-matrix eigenstates are 
retained. 

The DDMRG method and GPT work well together 
because of the following reasons. (1) It is not necessary 
to repeat DDMRG calculations because GPT does not 


impose self-consistency. (2) Data under open boundary 
conditions, which are obtained accurately in the DDMRG 
method, are used in GPT. The combined method can 
be regarded either as CPT where the DDMRG method 
is utilized as a cluster Green function solver or as the 
DDMRG method where momenta are interpolated based 
on GPT. Note that real-space cluster Green functions 
are used in the combined method [Fig. HJe)] in contrast 
with the random-phase approximation (RPA) from the 
decoupled-chain limit^ using DDMRG results [Figs. [2][c) 
andllKf)]— In this paper, the RPA from the decoupled- 
chain limit, which corresponds to the perturbation theory 
up to the first order with respect to interchain hopping, is 
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only used to explain how the spectral weights are shifted 
by interchain hopping from the decoupled-chain limit and 
to trace the origins of the dominant modes of 2D systems 
back to those of one-dimensional (ID) systems (Secs. 

Hirai and llirDll . 


C. Cluster-size effects 

As shown in Figs. HKc), [IKe), HKi), HKm), andUKn), the 
results obtained using (6 x 6)-site CPT are almost the 
same as those using (4 x 4)-site CPT. This implies that 
the overall spectral features will not change significantly 
if the cluster size is increased. 


III. SPECTRAL PROPERTIES 
A. Dominant modes 

The results for the 2D t-J model obtained using CPT 
are shown in Fig. [T] The overall spectral features can 
be explained in terms of the dominant modes [Fig. (He), 
modes i'-vi'] corresponding to those of the 2D Hubbard 
model [Fig. HKg), modes i-vi];^ whose origins can be 
traced back to those of the ID models [Fig. HKd), modes 
Fig. UKa), modes 1-5}^“— by considering how 
the spectral weights are shifted by interchain hopping 
from the decoupled-chain limit [Figs. |2l[c) and|2l[f); Sec. 
IIIID| i "^°d^i58i59 dispersing mode around (tTjTt) for 
a; > 0, mode i' (~ mode i), originates from mode 1' (~ 
mode 1, upper edge of the spinon-antiholon continuum). 
The mode around (0,0) for a; < 0, mode ii' (~ mode 
ii), primarily originates from mode 2' (ai mode 2, spinon 
mode), and the mode for slightly lower w, mode iii' (~ 
mode hi), primarily originates from mode 3' (~ mode 
3, holon mode). The mode spreading over the Brillouin 
zone in the large negative uj regime, mode iv' (« mode 
iv), primarily originates from mode 4' (Ri mode 4, holon 
mode called the shadow band). The mode bending back 
around (7r/2,7r/2) for w > 0 [r: upper edge of the band 
around (7r/2,7r/2) at half-filling], mode v' (r mode v), 
originates from mode 5' (r mode 5, antiholon mode). 
The flat mode, mode vi' (r mode vi), is dominant around 
(7r,0). 


B. Positive to 

Mode i' corresponds to the doping-induced states 
observed in high-Tc cuprates and in theoretical 
calculation a^dd^d^, 20-24,32,40,42,60- _76 controversial 

interpretations. The CPT results indicate that the 
energy of this mode at (tt, tt) [e(7r, tt)] does not reach 
zero even in the small-doping limit [Fig. Iljn)]. The 
extrapolated value of e(7r,7r) to J —>■ 0 [e( 7 r, 7 r)] behaves 
essentially as v^Us in the small-J/t regime, where Vs 



FIG. 2: (a) A{k,uj)t for U/t — 8 at 5 = 1/15 (r 0.07) 
for the lower Hubbard band (LHB) of the ID Hubbard 
model^ on a 120-site chain obtained using the non-Abelian 
DDMRG method with 120 density-matrix eigenstates. The 
blue numbers indicate dominant modes, (b) n{k) for (a) for 
tj < 0 [n_(A:)] (solid blue curve with hatches), w > 0 in the 
LHB (dotted red curve), the upper Hubbard band (dashed 
pink curve), and all oj [nt(fc)] (solid brown line).— (c) RPA 
results for the interchain hopping integral t± = t obtained 
using the data in (a).— (d)-(f) The same as (a)-(c) but for 
the t-J model at J/t = 0.5. (g) The same as the left panels 
of Fig. (He) but for the LHB of the 2D Hubbard model 
But U/t = 8.— The blue numbers indicate dominant modes, 
(h) n{k) for (g) with the same line types as in (b).— The 
solid green lines indicate u) = 0. Gaussian broadening with a 
standard deviation of O.lt is used. 


denotes the spin-wave velocity of the 2D Heisenberg 
model {vs R 1.18-\/2J [s^) [Fig. dKo)], as in the 2D 
Hubbard model.— In addition, the spectral weight for 
a; > 0 (FF+) behaves exactly as S [Fig. |Tl[i)]i^i2^ 
These results imply that mode i' continuously leads to 
the magnetic excitation of the Mott insulator, while its 
spectral weight gradually disappears toward the Mott 
transition. This feature is essentially the same as that in 
the ID and 2D Hubbard model s'^^d^i^^ and is consistent 
with a general argument in the small-doping limit. 
Thus, this feature will be related to transformation to the 
Mott insulator, which has low-energy spin excitation but 
no low-energy charge excitation (spin-charge separation), 
rather than being related to double occupancy. 

Because mode v' bends back around (7r/2, tt/2) for w > 
0 [Figs. [Ti;b),[Ti;c), and [TJe)], there are small intensities 
at momenta slightly away from (7r/2,7r/2) toward (tt, tt) 
as well as significant intensities around (7r/2,7r/2) for 
w R 0 near the Mott transition [Figs. [TJj) and[T][k)] — 
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The region surrounded by these intensities might be 
effectively regarded as a hole pocket. Regarding hole 
pockets, the momentum distribution function n_(fc) has 
been investigated ) a high momentum resolution 

is desired near the Mott transition at zero temperature. 
The CPT results indicate that n_ (fc) exhibits a small dip 
near (7r/2,7r/2) [Figs. [Iff) and[T][g)], for which mode v' 
would be more or less responsible4i‘^‘2& 


C. Small negative cj 

The bandwidth of mode ii' at half-filling [£( 0 . 0 )] 
[Fig. [Ha)] has been studied in relation to the hole 
motion in an antiferromagnet, and its J dependence 
(particularly the power-law behavior in the small-J/t 
regime) has been discussed Figure [TJo) 

indicates that it behaves essentially as \/2vs in the small- 
J/t regime, as in the 2D Hubbard model4Sii2iiS^ The 
bifurcation between this mode and mode iii' as well as the 
reduction in spectral weight just below the bifurcation 
point [Figs. IHb), (He), and He)] can be identified^ 
as the giant kink and waterfall behavior observed in 
high-Tc cuprates^^— for which various interpretations 
have been proposed^ ~^^i^°'^^d°'^^i’^^i^^~ -l^ Modes ii' and 
iii' correspond to spinon-like and holon-like branches 
observed in high-Tc cuprates^i respectively. 

The properties around (tt, 0) are primarily 
characterized by mode vi', which corresponds to the flat 
band observed in high-Tc cuprates and in theoretical 

calculations J^ d^d3i40,63-66,7i ,ra. 9 ^- 97 , 102 , 103 ^ Because this 

mode carries large spectral weights and exhibits an 
almost flat dispersion relation over a wide momentum 
range around (tt, 0), it significantly contributes to the 
main peak of the single-particle density of states p{(jj) 
[Figs. H^)“Hc)]- For instance, as shown in Fig. Hp)) 
the energy difference of the main peak of pioj) from the 
top of the band at half-filling [sp] [Fig. Hs-)] is primarily 
determined by that of this mode at (tt, 0) [£( 7 r,o)] [Fig- 
Ha)]- This energy difference, which can be regarded as 
a pseudogap in the small-doping limit, will be related 
to the antiferromagnetic fluctuation because it is almost 
proportional to J in the small-J/t regime [Fig. Hp)]) as 
in the 2D Hubbard model4aiS3,65 

Near the Mott transition, because mode vi' is located 
below w = 0, there is no mode crossing w = 0 along 
(0,0)-(7r, 0), and the spectral weights along (tt, 0)-(7r, tt) 
for w Ri 0 are also reduced significantly [Figs. Hb),Hc), 
He), Hj), and Hk)]*^ Thus, the spectral weights for 
a; Ri 0 essentially remain only around (7r/2,7r/2), which 
corresponds to the Fermi arc behavior observed in high- 

Tc cupratesi^i^ 

In the large-doping regime, the behavior becomes 
similar to that of noninteracting electrons after the 
pseudogap closes [the peak of p[ijj) or the flat mode 
crosses w = 0] [Figs. Hd),HF), andHl)]"^*^” 


D. Large negative ui 

Although mode iv' corresponds to mode iv, its spectral 
weights around (0,0) and (tt, tt) are significantly smaller 
and larger, respectively, and its lo values around (0,0) 
and (tt, tt) are significantly higher than those of mode 
iv near the Mott transition [Figs. Ht>), H^), H®), 

Hg)] - Similar features have also been obtained primarily 
at half-filling using exact diagonalizationiSi and recently 
using CPT independently!^ Here, we interpret these 
features as a consequence of the restriction on spectral 
weights. 

For the electron operators with the no-double- 
occupancy constraint, the spectral weight for each 
k [nt(fc)] is reduced to (l-|-i5)/2 [Figs. Hf)“Hh) and 
He)]- As a result, the spectral weight for w < 0 [n_(fc)] 
around fe = 0 in the t-J model becomes significantly 
smaller than that of the Hubbard model near the Mott 
transition [Figs. Hs) and H^); Figs. He) and H^)]- 
Because the spectral weight for w < 0 (IF-) is equal to 
(1 — 6)/2 in the t-J and Hubbard models, n-{k) around 
fc = TT in the t-J model becomes larger to compensate for 
the reduction around fc = 0. Here, 0 and tt respectively 
indicate 0 and tt for ID systems and (0, 0) and (tt, tt) for 
2D systems. 

In the RPA from the decoupled-chain limit 
[G“^(fc,a;) = G^^{kxTU}) — t_L(fc)]^ the spectral weights 
shift upward and downward in the momentum regime 
for t^{k) > 0 and t_L(fc) < 0, respectively ; 
where Gi-D{kx,oj) and t±(k)(= —2t±cosky) denote the 
Green function of a chain and Fourier transform of the 
interchain hopping integral, respectively. In addition, 
the mode shift becomes large if the mode carries large 
spectral weights in the large-|ti(fc)| regime 
Because the spectral weight around fc = 0 for a; < 0 in 
the ID t-J model is significantly smaller, the downward 
spectral-weight shift around (0, 0) is smaller [Fig. HO] 
than that of the Hubbard model [Fig. H^)]- 

This argument explains why the uj value around (0,0) 
of mode iv' is higher than that of mode iv near the Mott 
transition [Figs. H^) andHg)]- For (7r,7r), a similar 
argument can also explain that the u value around (tt, tt) 
of mode iv' is higher than that of mode iv. Near the 
Mott transition, however, the w values of modes 4, 4', 
iv, and iv' around fc = tt appear to be better explained 
by considering that they are almost the same as those 
around fc = 0 [Figs. H3')“Hc),H6),H9'),Hd), andHg)]- 
At half-filling, the oj value of mode 4 at fc = tt is exactly 
the same as that at fc = 0 because the dispersion relation 
of the holon mode becomes symmetric with respect to 
fc = 

For the differences in the large negative uj regime, the 
neglect of the three-site term is relevant because 

this term contributes to hopping. 
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IV. DISCUSSION AND SUMMARY 

In this study, the single-particle spectral properties of 
the 2D t-J model for J/t r; 0.5 near the Mott transition 
were investigated. In contrast with conventional exact 
diagonalization studies, where the spectral weights were 
calculated only at available k and S points depending on 
the cluster size and boundary conditions, the spectral- 
weight distributions for continuous k and ui near the 
Mott transition with small intervals of 6 were calculated 
using CPT, and how the peaks of the spectral weights 
form dominant modes and how the modes transform 
to those of the Mott insulator as 6 gradually decreases 
were clarified in the 2D t-J model. In addition, through 
comparisons of the (4 x 4)-site CPT results with the 
(6 X 6)-site CPT results obtained by combining the non- 
Abelian DDMRG method, the cluster-size effects were 
confirmed to be small enough to allow discussion on the 
overall spectral features. 

Furthermore, the natures of the dominant modes in 
the 2D t-J model were clarified by investigating the S 
and J/t dependences of the characteristic energies and 
spectral weights and by tracing the origins of the modes 
back to those of the ID models. In terms of the dominant 
modes, various anomalous spectral features observed in 
high-Tc cuprates, such as the doping-induced states, flat 
band, pseudogap, Fermi arc, spinon-like and holon-like 
branches, giant kink, and waterfall behavior^^— were 
collectively explained in the 2D t-J model, as in the 2D 
Hubbard model.— 

The results for the natures of the dominant modes 
imply that these spectral features are primarily related 


to the proximity of the antiferromagnetic Mott insulator, 
which has a low-energy spin-wave mode^SS but no low- 
energy charge excitation, and to the existence of different 
energy scales that characterize the bandwidths of the 
dominant modes rather than to double occupancy, which 
is completely removed in the 2D t-J model. 

Although the spectral features of the 2D t-J model 
in the small-|a;| regime are similar to those of the 2D 
Hubbard model, significant differences appear in the 
large negative lu regime around (0, 0) and (tt, tt) near the 
Mott transition for the modes which can be considered 
to primarily originate from the ID holon modes. In this 
study, the differences were interpreted as a consequence 
of the restriction on spectral weights. 

Because of the limited resolution, this study could not 
clarify the properties in the very small-|a;| regime, such 
as the nature of the excitation in the small-|a;| limit, the 
accurate gapless points of the single-particle spectrum, 
and the presence or absence of a long-range order in the 
ground state, although a superconducting ground state 
has been suggested in a considerable number of studies 
for the 2D t-J model— Further studies are 
needed to clarify more details and how the anomalous 
features are related to high-Tc superconductivity. 
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